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Abstract 

Atom and nanoparticle arrays trapped in optical lattices are shown to be capable of sustaining 
collective oscillations of frequency proportional to the strength of the external light field. The 
spectrum of these oscillations determines the mechanical stability of the arrays. This phenomenon 
is studied for dimers, strings, and two-dimensional planar arrays. Laterally confined particles free 
to move along an optical channel are also considered as an example of collective motion in partially- 
confined systems. The fundamental concepts of dynamical response in optical matter introduced 
here constitute the basis for potential applications to quantum information technology and signal 
processing. Experimental realizations of these systems are proposed. 
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INTRODUCTION 



Periodic standing lierht patterns created by coherent laser beams offer tfie possibility to 
trap objects ranging from single atoms [jj to micron-size particles |2|, |3[, forming ordered 
arrays that can be viewed as artificial crystals with the periodicity of the optical field. 
Starting from the equilibrium configuration of an ordered array with one particle per lattice 
site, the perturbation produced by slightly displacing one of the particles will be transmitted 
to the others via their mutual electromagnetic interaction, which is a complex function of the 
particles positions. This results in collective motion of the particles that resembles phonon 
vibrations in crystals and that we regard as collective oscillations in optical matter. 

In this context, the interaction of ultracold atoms with optical lattices has recently at- 
tracted considerable attention because of its potential use for quantum information tech- 
uology. ,n which the reahzafo. of trapp.ug of no more than one atom per latfee site has 
been an important benchmark [l)]. Using light tuned close to some atomic resonance, strong 
light-induced interaction between atoms trapped in different sites can take place even for 
commonly employed lattice periods of the order of 1 /im, as suggested by recent observa- 



a. 



These interatomic 



tions in optical cavities under atom-cavity resonance conditions 
interactions are the driving force of the collective oscillations considered here, which for 
reasonably-cold atoms offer a potentially practical realization of quantum gates 6|] involving 
large numbers of qubits. 

Early attempts to bind small particles using light forces led to the 

development of optical tweezers , capable of trapping and aligning objects ranging from 



micro-organisms 



10| to metallic nanoparticles ll|, ll2| . Furthermore, holographic tweezers 



have been developed to trap tailored arrays of micron-size objects [13|, that offer an 
excellent playground to test many of the concepts discussed below. Manipulation of micro- 
particles using plasmons has been recently demonstrated as well yL5i], whereas fine tuning of 
nanoparticle positions has been theoretically proved to be realizable by coupling to plasmonic 
nanostructures 16|. 

Mutual inter-particle interaction induced by external illumination has been shown to 
lead to interesting effects such as configurational bistability in pairs of spheres 17| and 
crystallization in linear particle arrays trapped at the focus of counter-propagating lasers 
181] . A close relative of the oscillations studied here are those of colloidal crystals 19| . 
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mediated by electrostatic interaction similar to plionons in ionic crystals. 

Here, we examine the collective oscillation modes of atoms and particles trapped in optical 
lattices. For atomic lattices, strong inter-atomic interaction effects are predicted in the 
oscillation spectrum near atomic absorption resonances. The particle dimer is studied first 
as a tutorial example. For extended systems, long-range 1/r dynamical interaction between 
polarized atoms or particles is shown to give rise to soft modes observed in particular in 
one-dimensional (ID) periodic arrangements. Finally, periodic arrays confined within a ID 
well are shown to exhibit imaginary collective-motion eigenfrequencies for some ranges of 
the particles spacing, thus revealing instabilities that preclude structural transformations. 



THEORETICAL BACKGROUND 



We start by analyzing collective motion of small particles such as atoms, molecules or 
nanoparticles that respond to electromagnetic fields basically as induced dipoles driven 
by their polarizability a{uj), where u is the light frequency. The time- averaged force 
acting on one of such particles in the presence of an external electric field W^^{r,t) = 
2 3?{E'^^*(r) exp(— icjt)} is readily obtained from the integral of Maxwell's stress tensor on a 
small sphere surrounding the particle. One finds 
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where I labels Cartesian coordinates, ^(r) = — 3?{Q;}|E'^^*(r)| acts as an effective potential 
(i.e., the particle can be trapped in regions of low or high electric field strength when is 
negative or positive, respectively), the second term in the right hand side of Eq. ([1]) describes 
the so-called radiation pressure, and we use Gaussian units throughout this paper. 

Focusing on atoms, the frequency-dependent atomic polarizability near a resonance at 
frequency ujq with no decay channels other than radiative takes a simple form compatible 
with the optical theorem condition Q{ — l/a} = 2k^/3 2l| : 

^ 3c3r/4^3 



a{uj) 



3c3r/2cj2 



(2) 



Uq — u'^ — iFtu^/tUo uq — u — ir/2' 
where k = u/cis the momentum of light in free space, F is the resonance frequency width, 
and the last approximation is valid for T <^ ujq and u near the resonance. Furthermore, our 
analysis assumes low-enough light intensities to exclude saturation effects near resonance, 
which would require to go beyond linear response approximation. 



When several particles placed in vacuum are considered, their induced dipoles express 
the response to the external field plus the field scattered by the other particles, that is. 



a 



E^^*(r,) + EG'(r,-r,Op,' 



(3) 



where tj labels particle positions and an exp(— itut) time dependence is understood. The 
3x3 matrix 



G{r) = (A;2 + V ® V) = A 



B 



(4) 



with coefficients 



A 
B 



{krY + ikr — 1 



(kr) 



3ikr — 3 



(5) 
(6) 



describes the electric field produced at the position r by a dipole at the origin. All particles 
are considered to be equal for simplicity. The force acting on the particle at rj is then 
obtained from Eq. ([T]) by adding the scattered electric field J2j'ytjG{r — rj/)pj' to E'^''*(r). 



Up to here, our analysis follows previous developments o: 
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optical forces acting on small 



23| . This type of formalism can 



particles described through their dipolar response 
only be applied to particles that are small with respect to both the wavelength and their 
separation. However, its extension to include multipolar terms in the interaction between 
neighboring particles has been successfully carried out both for axially-symmetric particles 



using 3D multipo 



multipoles 



26|. 



ar expansions 2J] and for objects with translational invariance using 2D 



It is easy to see that the dynamics can no longer be obtained from an effective potential 
V, even in the absence of radiation pressure. This is a manifestation of the fact that we are 
dealing with open systems in which the photon bath provided by the external field can add 
(remove) energy to (from) the motion of the particles. 

OPTICAL DIMER 



We shall discuss first a simple atomic dimer system that illustrates the effect of strong 
inter- atomic interaction near an absorption resonance. Optical biding in a dimer has been 
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extensively studied both theoretically [27| and experimentally [2I, Q, but here we shall focus 
on the dynamical aspects of such system. 

An optical lattice will be considered to be formed by three pairs of equal-intensity counter- 
propagating lasers that define an orthogonal xyz frame, with their linear polarizations as 
shown in the lower inset of Fig. [TJ The external electric field produced by the lasers is then 
given by 



E' 



ext 



>, t) = 4 3ft e"'"^*} [?,m{kz) x + sin(A;a;) y + ?,in{ky) z] 



where Eq is the electric field amplitude of each laser. For uj > uoq [see Eq. ([2])], the atomic 
polarizability is essentially negative, as shown in the upper inset of Fig. [H which represents 
the polarizability of Eq. ([2]) scaled to the wavelength (a/A'^), so that the two atoms under 
consideration can be trapped at contiguous minima of the electric field (light regions in the 
lower inset). 

The oscillation frequencies, as derived from Newton's equation for infinitesimal displace- 
ments from equilibrium positions, reduce after lengthy, straightforward algebra from Eqs. 
m and (BD to 



± 



\E.\k 



3ft 



±A 



(7) 



and 



\Eo\k 



M \ 



3ft 



-1 



a 



-1 



±iA-B) 



(8) 



where M is the atomic mass, A = — 8(7r^ + ivr — 1)/A^, and A — B = 16(i7r — 1)/A^, according 
to the definitions of Eqs. ([5]) and ([6]) for a dimer spacing r = A/2. The light-frequency 
dependence of these oscillations is shown in Fig. [H along with schematic representations of 
directions of motion for the modes, which illustrate how different orthogonal directions are 
decoupled. Modes of frequency ilf and ilf and doubly and singly degenerate, respectively, 
and they are represented by solid and dashed curves. 

The dimer oscillation frequencies deviate considerably from the single atom case (dotted 
curve) as a result of inter-atomic interaction driven by the ujq resonance. Interestingly, 
this interaction pushes the range of instability of the dimer (signaled by some imaginary 
eigenfrequency Q) slightly towards the uj > ujq region as compared to the single atom. 
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FIG. 1: Oscillation frequencies of two identical atoms trapped in contiguous wells of an optical 

lattice as a function of light frequency uj near an atomic resonance ujq. The lattice is set up by 

three pairs of counter-propagating lasers of direction and (linear) polarization as shown in the 

lower inset. The wells are a distance A/2 apart, where A = 2ttc/u) is the light wavelength. The 

atomic polar izability a (upper inset) is obtained by assuming that radiative decay is the dominant 

contribution to the resonance width T [see Eq. ([2])]. The oscillation frequency Vt is normalized 

using the atom mass M and the intensity of each beam /. The oscillations take place around 

equilibrium positions corresponding to the vanishing of the electric field strength (yellow regions 

in lower inset). Six oscillation modes are obtained: two pairs of doubly-degenerate modes (solid 

curves) and two nondegenerate modes (dashed curves). The relative directions of motion of the 

atoms in the dimer are indicated by arrows for each of the modes. The oscillation frequency of a 

single trapped atom is given for reference (dotted curve). 
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Optical forces scale linearly with light intensity, and therefore, the oscillation frequencies 
are proportional to the strength of the applied electric field. This allows us to represent those 
frequencies normalized to the light intensity flux per laser, / = \Eq\'^c/2ti. For instance, 
single Rb atoms (M = 85.46 amu) near their 780 nm resonance (F = 27r x 4 MHz 28|, |29| |) 
will oscillate with frequency Vt ^ 0.36 MHz for / = 10 mW/cm^ lasers tuned to u = ti;o + 2F, 
which corresponds to a temperature of ~ 17 fiK. Notice that f2 <^ F, so that Doppler shifts 
can be safely neglected at this relatively high laser intensity, for which the atomic motion 
can be described classically, although saturation effects could be an issue Q]. 



PERIODIC ARRAYS 



Oscillations in periodic configurations of extended arrays can be studied by introducing 
small perturbations in Eq. for specific momentum q, corresponding to a displacement of 
every particle j around its equilibrium position given by Uj = Uqexp(iq ■ rj). For small 
perturbations, the electromagnetic forces scale linearly with the displacement Uq, so that 
one obtains an equation of motion of the form 

K{Sq}Uq = -Mfi^Uq, (9) 

where M is the particle mass. Three branches of the oscillation frequency fiq are obtained 
in general out of the eigenvalues of the force constants matrix Sq. 



ID arrays 

We shall first focus on blue-detuned systems, in which the particle polarizability has 
negative real part, so that the unperturbed particles sit at light intensity minima. Then, 
the force constants matrix reduces to 

Sq = \Eo\'e C+^^C, 

a H 

where 

G'q = Ee-^"-^^G'(r,) (10) 
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is the sum of the dipole-dipole interaction extended over lattice sites r^, and C is a di- 
mensionless matrix that depends on the specific orientations and polarizations of the lasers 
setting up the optical lattice. 

In particular, we shall consider a ID infinite string formed by the same kind of blue- 
detuned atoms as in Fig. [H with one atom trapped at each minimum of the electric-field 
along the (111) direction of the same optical lattice. The oscillation frequencies are then 
labeled by the longitudinal momentum q. Displacements parallel (||) and perpendicular (±) 
to the array turn out to be decoupled. We find the following analytical expression for the 
oscillation frequencies, similar in structure to Eqs. ([7]) and ([H]): 



\Ea\k 



3? 



-1 



where a =||, -L in each case, the lattice sums reduce to 



= 4^cos(gaj)^— y^(l - iA;aj) (11) 
and 

Gf = 2Y,cos{qaj)—-^[{kajf + ikaj - 1], (12) 
j=i W) 

and a = \/3X/2 is the period of the array, so that ka = y/Sn. 

Fig. [2] shows results obtained using this formalism within the first Brillouin zone in q 
space. Soft modes of zero frequency can be observed in the degenerate oscillations along 
directions transversal to the string, giving rise to large group velocities. This results from 
phase accumulation in the long-range inter-atomic interaction within Eq. f|T2|) . which di- 
verges as In ± ga| near =pga = = (2 — \^)tt. This logarithmic divergence is obviously 
destroyed when the string is finite, regardless its length, and therefore, it cannot lead to 
permanent distortions of the lattice from the equidistant array configuration. 



Beyond the dipole approximation in 2D arrays 

For particles of non-negligible size compared to the wavelength, multipoles beyond dipoles 
may become important. This situation can be described by a relation similar to Eq. Q, 
where p is then understood as a vector of multipolar amplitudes and a is the multipolar 
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FIG. 2: Dispersion relation of the oscillation modes in an infinite periodic string of atoms located 
at the minima of electric field strength along the (111) direction under the same conditions of 
illumination as in Fig. [1] (see lower inset in that Fig.; the coordinate axes are taken parallel to 
the laser beams). The nearest-neighbor distance is a = \/3A/2. The light frequency is tuned 
to u = + r? using the same notation and atomic polarizability as in Fig. [TJ The oscillation 
frequency is given as a function of momentum q along the string within the one-dimensional first 
Brillouin zone. Oscillations parallel to the (111) direction (broken curve) are decoupled from those 
along perpendicular directions (solid curve). 



scattering matrix [3(|. Analyses along these lines have been previously offered by the au- 
thor for momentum transfer from a fast electron to a particle 311] and for neighboring effects 



in the interaction force between two arbitrarily-shaped particles [2J]. Higher-order multi- 
poles beyond the dipole yield a more complex expression for the force, but the quadratic 
dependence on the field exhibited by Eq. ([1]) is still maintained. 

These conditions have been considered in Fig. [3] for a two-dimensional (2D) square array 
of silica beads (see upper inset) trapped in a field given by 

E""*(r,t) = 8 3ft{^o e"^"*} { [sin(Qx) x + sin(gy) y] sin (A;, 2) cos^ 

+ [cos(Qx) + cos((5y)] z cos(A;2-2) sin 6^1, 

where Q is the angle between the beam directions and the normal to the plane of the array, 
Q = ksinO, and kz = kcosO. This system is experimentally challenging because it requires 
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FIG. 3: Dispersion relation of oscillation modes in an optical lattice of silica spherical particles 
(e = 2.1) in vacuum. The lattice is set up by four pairs of counter-propagating linearly-polarized 
beams with incident magnetic field in the plane of the particles. The inset shows the particles in 
their equilibrium positions forming a planar square lattice of nearest-neighbor distance X/s'm9, 
where 9 = 50° is the angle formed by the light beams and the normal to the plane of the array. 
Continuous, dashed, and dashed-dotted thick curves correspond to the three oscillation modes for 
each value of the parallel momentum along the excursion rXMT within the first Brillouin zone. 
The spheres diameter is 240 nm, the light wavelength is A = 550 nm, and the flux of each laser 
is 100 W/cm^. These results are obtained with inclusion of higher-order multipoles beyond the 
dipole. Calculations in the dipole approximation (thin dotted curves) are shown for reference, with 
the polarizability obtained from the electric dipole Mie coefficient as a = (3/2 A;'^) tf 3C|. 



to trap particles in vacuum, for which pioneering results were reported three decades ago by 
Ashkin and Dziedzic 3^, who achieved optical levitation in high vacuum. 



Three oscillation frequency bands are obtained and represented for an excursion along 
symmetry points within the first Brillouin zone of the reciprocal lattice in momentum space. 
One of the bands corresponds to nearly-independent-particle motion perpendicular to the 
plane of the array and is almost dispersionless (dashed-dotted thick curve). The remaining 
two bands (solid and dashed thick curves) are degenerate at the F and M points, which is 
a condition imposed by symmetry. Unlike the previous examples, the dielectric particles 
of Fig. [3] are attracted towards maxima of the light intensity (this is the equivalent of red- 
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detuned atomic resonances). The effect of multipolar interaction is quantitatively significant 
(cf. oscillation modes calculated in the dipole approximation, represented by thin dotted 
curves in Fig. [21 incidentally, these latter dipole calculations should describe well red-detuned 
atoms trapped in this optical-lattice geometry). 

Low-frequency modes arising from long-range interactions are visible in Fig. [3] as pro- 
nounced dips that become very sensitive to finite array boundaries and size distribution of 
the particles. These effects have been phenomenologically described in Fig. [3] through a small 
attenuation in the interaction between distant particles represented by a dielectric function 
equal to 1 + 0.03i in the surrounding medium, which gives an interaction decay-length of 
~ 8 times the nearest-neighbor distance. 

PHASE TRANSITIONS IN ID OPTICAL MATTER 

The above analysis has been limited to particles trapped at the sites of 3D lattices and 
capable of oscillating around their equilibrium positions. We shall now explore a different 
system comprised by particles that are confined along two spatial directions but free to 
move along the remaining third direction. This scenario is presented for instance in particles 
trapped inside a cylindrical optical cavity (e.g., a hole of a photonic crystal fiber), and also 
in particles confined to a ID well of a 2D optical lattice like that of Fig. Il](a), in which the 
external field is 



where is the field amplitude of each laser. The particles will be assumed to have positive 
polarizability a (i.e., the confinement occurs in regions of maximum electric field) and to 
be periodically spaced. This is a configuration of equilibrium for an infinite chain, in which 
the force acting on the particles vanishes. We shall then analyze the spectrum of collective 
oscillations within the first Brillouin zone of the reciprocal lattice and explore the stability 
of that equilibrium. 

The self-consistent electric field amplitude acting on a given particle reduces from Eq. 
([3]) to the analytical expression 




[cos{kz) + cos(%)] X, 
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FIG. 4: Oscillation modes and stability of a linear array of particles trapped in a two-dimensional 
optical lattice set up by two pairs of counter-propagating beams with polarizations as shown in 
(a). The array is assumed to be periodic, infinitely long, and surrounded by vacuum. A typical 
spectrum of oscillations is presented in panel (b) for specific values of the free-space momentum 
of the trapping light k and the particle polar izability a, normalized to the lattice spacing a (see 
text insets). The spectrum is given within the first Brillouin zone of momentum along the array 
q, with normalization of frequency and momentum as in Fig. [5J The contour plot of panel (c) 
represent, squared oscrOation trequeneres as a tunct.on of , and t tor „ = l/(<.„-' - 2itV3). where 
ao = O.OOla^ is the electrostatic polarizability (this prescription for a is consistent with the optical 



where Gg=o [see Eq. ffTTl) ] is the sum of the fields scattered by the rest of the particles, and 
only the exp(— icut) component is considered. Small displacements out of equilibrium along 
small vectors Uj, as defined at the beginning of Sec. , will then modify the self-consistent 
field in the vicinity of every particle j to 

from which we discount the field scattered by that particle. The field perturbations can be 
conveniently obtained by expanding G{r) [see Eq. (jl])] in Taylor series and by working in 
Fourier space q according to the definition 



-iqaj 



similar to Eq. ( fTOl) . We find 



9qE 
~ — Gn 



U" 



and 



where || and _L indicate directions parallel (x) and perpendicular (y, z) to the array, respec- 
tively, Gj and G^ were defined in Eqs. (fTTl) and (fT2l) . and 

oo gifcai 

9q = 4i5Isin(gaj)^—y^[(A;aj)^ + 3ikaj - 3]. 

Finally, the oscillation frequencies flq are obtained from Newton's equation, similar to Eq. 
(Q, but with a more complex expression for Sq. Like in the system of Fig. [21 oscillations 
along parallel and perpendicular directions with respect to the array are decoupled and their 
frequencies reduce to 



2\aE\ 
M 



3? 



9a 



1/a — Gq 



1/2 



(13) 



and 



2\aE\ 
M 



3ft 



1/a - G^ 



M 



1/2 



(14) 
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respectively, where 



oo 



,ikaj 



Hi =4 5][cos(gaj) - 1] 



[i{ka]f - 5{kajf - Uikaj + 12] 



i=i 



and 



oo 



,ikaj 



H^ = 8Y,[cos{qaj)-l] 



[{kajY + 3iA;aj — 3]. 



Perpendicular modes span a nearly-flat band around a high central frequency determined 
by the confining transversal potential [see the right-most term inside the square brackets 
of Eq. (IT^ ]. The stability of the array will be exclusively determined by the band of 
longitudinal modes. A typical dispersion relation is shown in Fig. IH^b) for specific values 
of the polarizability (see inset) and spacing-to-wavelength ratio (a/A = 1.1). In contrast 
to the system of Fig. [21 the longitudinal band is acoustic, so that the oscillation frequency 
vanishes in the g ^ limit (notice that = = {] in this limit), standing for a rigid slow 
translation of all particles along the ID trapping well. 

The stability of this type of array is analyzed in Fig. Hl^c), which represents the oscillation 
spectrum given by Eq. (IT^ as a function of lattice spacing. White regions correspond to 
modes of imaginary frequency that describe motion beyond the validity of the harmonic 
approximation. The presence of these modes signals structural instabilities for some ranges 
of the period of the array, which must rearrange itself by creating defects with different 
particle separation. However, there are certain ranges of spacings that are stable, centered 
around the distance of equilibrium for a particle dimer in the channel (see horizontal arrows). 
This analysis of the array stability is also valid in the presence of friction forces, like those 
that would show up if the particles were immersed in a fluid. A practical realization of such 
arrays using holographic tweezers [is], [3] combined with transversal confining light would 
allow one to study their stability and to assess mechanical properties like the compressibility 
by fixing the position of the particles at the end of a finite array. 

CONCLUSION 

Atoms and nanoparticles trapped in optical lattices have been shown to exhibit collective 
oscillations, the spectra of which reveal complex patterns that can be controlled by exter- 
nal illumination conditions. The oscillation frequency increases with the amplitude of the 
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binding lasers, thus adding an extra degree of freedom that allows obtaining frequencies in 
the kHz-MHz range in the case of atomic lattices. These oscillations constitute a genuine 
form of dynamical collective behavior in optical matter encompassing complex phenomena 
such as soft-modes and lattice phase transitions. Furthermore, our analysis of the oscillation 
spectra has been shown to resolve the question of the stability of the arrays. 

The results presented here could be useful to design the following specific experiments: 

• Collective motion in optical lattices of trapped atoms, which should be observable in 
their response to low-frequency radiation within the noted oscillation frequency range. 

• Stability of ID nanoparticle chains trapped along one channel of a 2D optical lattice. 
Initial particle trapping is possible using holographic illumination p^], superimposed 
to the optical lattice, and the stability could be proved by varying the relative intensity 
of holographic and optical-lattice laser beams. 

• Phase transitions in optical matter, using a similar setup as in the previous point, but 
holding the ends of a chain at specific locations. Transitions should occur as these 
ends are displaced. 

Finally, there are a number of open questions posed by these novel oscillations that will 
require separate developments, including the extension of this work to lower laser intensi- 
ties and to conditions for which the oscillations behave as quasi-particles in atomic lattices. 
Interesting phenomena could also result from saturation effects at large laser intensities, 
whereby the atomic response is far from the linear regime. Equally important, thermody- 
namics of open systems like the arrays of Fig. H] is expected to yield new properties associated 
to physical quantities such as the compressibility when exerting pressure on finite strings 
or the specific heat associated to the oscillations. These are challenges that will ultimately 
determine the applicability of collective behavior of optical matter to fields such as quantum 
information technology and signal processing. 
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